Spin transport in coupled spinor Bose gases 
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We report direct measurements of spin transport in a trapped, partially condensed spinor Bose gas. 
Detailed analyses of spin flux in this out-of-equilibrium quantum gas are performed by monitoring 
the flow of atoms in different hyperfine spin states. The main mechanisms for motion in this 
system are exchange scattering and potential energy inhomogeneity, which lead to spin waves in the 
normal component and domain formation in the condensate. We find a large discrepancy in domain 
formation time scales with those predicted by potential-driven formation, indicating strong coupling 
of the condensate to the normal component spin wave. 

PACS numbers: 03.75.Kk, 03.75.Mn, 51.10.+y, 75.30.Ds 



An atomic spinor Bose gas, comprised of multiple spin 
species, exhibits rich dynamics not present in systems 
without the added spin degree of freedom. Additionally, 
a partially condensed system, with a Bose-Einstein con- 
densate (BEC) in thermal equilibrium with an uncon- 
densed component, adds immensely to the complexity. 
Such a gas can exhibit effects commonly associated with 
bosonic systems, including superfluidity and macroscopic 
phase coherence, as well as effects more typically asso- 
ciated with fermionic systems, namely spin waves and 
domain formation. This Bose ferromagnet system is an 
interesting system in which to study out-of-equilibrium 
quantum effects resulting from a coupling between the su- 
perfluid and normal components of the gas, particularly 
in light of recent interest in the coexistence of superflu- 
idity and ferromagnetism in solid-state materials 
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Recent work with tunable interactions in lattice-based 
spinor systems demonstrated phase diffusion 2] and spin 
squeezing [3| , and there has been much theoretical inter- 
est, particularly in strongly interacting spinor systems. 
However, even so-called weakly interacting systems with 
static interactions have shown complicated and surpris- 
ing behavior, and these systems continue to be of interest. 
Much work has been done with spinor systems at the two 
limits of the temperature spectrum - in nearly pure con- 
densates [J, [fj] and in thermal gases [a, [D], but due to 
the complexity of the system there have been relatively 
few studies ofpartially condensed spinor systems at finite 
temperature |8l4l3| . Early experiments in finite tempera- 
ture systems demonstrated condensate melting driven by 
normal component decoherence [10(, superfluid- normal 
component phase lockin g |l l|. and thermal decoherence 
in coherent spin mixing [12j. More recently, remarkable 
spin textures were observed in an optically confined Bose 
gas at finite temperature [13|. In that work, and sev- 
eral others using pure condensates [g, [l2|, interconver- 
sion between spin states was one of the dominant effects 
in the dynamics. Although interconversion leads to rich 
dynamics, individual spin populations are not conserved, 



which limits the type of information that can be obtained 
about the non-equilibrium dynamics of those systems. 
This work differs in that spin populations are strictly 
conserved, allowing for direct observation of spin-state 
flow during the unfolding spin kinetics. 

In Ref. [ll| . condensate- normal component coupling 
was inferred from anomalously fast domain formation as 
well as phase uniformity across both components. Here 
we confirm these observations with quantitative spin 
transport measurements that definitively rule out the 
usual mechanism of condensate domain formation [J, Il4j . 
As atoms confined in a magnetic trap evolve, their spins 
are altered by interactions with other atoms and by the 
confining potential itself. By monitoring the spatiotem- 
poral flux of atoms in the individual spin states forming 
the spinor, we follow energy transport driven by the spin 
kinetics. In particular, we measure the spatial distribu- 
tions of the spin species and derive spatially dependent 
velocities for each spin state from the profiles, indicating 
a spin flux. The mechanism for generating this spin flux 
is examined by comparing the measured velocities with 
velocities predicted from inhomogeneous transverse spin 
profiles, showing a large discrepancy, with the spin wave 
accelerating the condensate to speeds nearing the speed 
of sound. This velocity discrepancy can be reconciled 
through the coupling of condensate domain formation to 
the normal component spin waves. The presence of the 
condensate can be interpreted as an increased spin stiff- 
ness in the coupled system, as described in [9(. 

The spinor system described herein is comprised of 
two hyperfine spin states of the 5Si/ 2 manifold of 87 Rb 
atoms. These two states, |1) = \F — l,mp = —1) and 
|2) = \F = 2, mp = 1), form a pseudospin doublet, which 
may be thought of merely as a spin-1/2 system. Bloch 
vector notation is used to describe the spinor [15j; the 
coordinate axes of the Bloch sphere are taken to be u, v 
and w. The "longitudinal" spin component S w is the 
population inversion, and the "transverse" spin compo- 
nents S u and S v are the coherences, with the angle in the 
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FIG. 1: Typical radially averaged density profiles at T/T c = 
0.8 for the spin states |1) (dashed line) and |2) (solid line) af- 
ter a spin wave oscillation and domain formation have begun. 
The |1) state forms domains surrounding the |2) state, which 
resides in the center of the condensate (40 /im Thomas-Fermi 
radius). The inset image is a sum of the |1) and |2) images 
(expanded radially), showing no measurable total density fluc- 
tuation despite the spin dynamics. 



u — v plane <fi being the relative phase between states |1) 
and 1 2) in a coherent superposition. The wave-function 
of this system can be written \tp) = (|1) + e ,; ^|2))/\/2. 

Spin kinetics in a relatively pure, trapped spinor con- 
densate are dominated by the formation of segregated 
magnetic spin domains [4J. Atom-atom interactions, 
characterized in the ultracold scattering limit by the s- 
wave scattering length, a, are responsible for domain for- 
mation. The criterion for in-trap domain formation is 
that a spinor condensate's interspecies scattering lengths 
for the spin states must be large, positive, and not equal 
to the intraspecies scattering lengths 14 1. Domain for- 
mation does not significantly alter the overall density dis- 
tribution of the condensate, merely the density distribu- 
tions of the individual spin states (Fig. Q] inset). 

The spin flux needed to form condensate domains is 
produced by spatial gradients in the relative phase of the 
spin-state superposition. A spatial gradient in the phase 
of the superposition is a higher-energy configuration, and 
the condensate will redistribute energy by "untwisting" 
itself. More precisely, a gradient in the condensate phase 
leads to a velocity, u, as follows 16]: 



u{f) = (ft/m)V0(r), 



(1) 



where 771 is the atomic mass. Gradients in 4>{r) arise from 
spatial inhomogeneities in the potential energy, primar- 
ily from inhomogeneities due to the density-dependent 
mean- field shift (see Fig. [2]). In the trapping geometry 
used here, the typical time scale for domain formation is 
~ 280 ms for a nearly pure condensate. 

The experimental setup used for this work consists of 
87 Rb confined in an Ioffe-Pritchard magnetic trap similar 
to that described in [l7J. The trap is cylindrically sym- 
metric, with trapping frequencies of {6.7, 236, 236} Hz. 
Because the magnetic trap is highly anisotropic, all of 
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FIG. 2: Radially averaged phase profiles at t w ms (dashed 
line) and t = 70 ms (solid line), both at temperature T/T c = 
0.8, obtained using Ramsey spectroscopy. Initially, the phase 
is uniform across the cloud. As time evolves, spatial inhomo- 
geneities develop, largely due to the mean-field shift, which is 
dependent on the local density. The large phase bump in the 
center of the cloud shows the presence of the dense BEC. 



the significant dynamics occur only in the loosely con- 
fined axial direction, z, and the density distribution may 
be radially averaged and effectively considered to be a 
one-dimensional distribution (see Fig. [I}. The tempera- 
ture (and thus the fraction of atoms in the condensate) is 
controlled with radio frequency evaporative cooling from 
the critical temperature, T c , down to T/T c < 0.3. The 
two spin states are connected by a two-photon microwave 
coupling |6j. The initial conditions for the experiments 
described herein are immediately following the creation 
of an equal coherent superposition of the |1) and |2) states 
via a 7r/2 pulse. The initial condensate size is kept fixed 
at a small number, Nq ~ 6 x 10 4 atoms, so as to minimize 
the effects of loss via dipolar relaxation in |2). 

The spin flux is characterized by measuring the veloc- 
ities of the two spin states. This measurement is accom- 
plished as follows. The initial conditions are realized by 
creating a partially condensed condensate in a superposi- 
tion of the two spin states. Following an evolution time t, 
the spatial distribution of the population in either of the 
spin states can be measured (using destructive absorp- 
tive imaging). By analyzing the time-rate of change of 
the density profiles of each spin state (Fig. [T]), a relative 
velocity profile of the spin components is found. This 
relative velocity, u(z,t) = U2(z,t) — U\{z, t), is derived 
from drii/dt by integrating the continuity equation, 



drii d 

at dz K ' 



(2) 



where rii(z, t) and Ui(z, t) are the radially averaged den- 
sity and velocity for state \i). The final term accounts 
for density-dependent losses. The inset of Fig. [3] shows 
a typical velocity profile. Because the atoms' velocities 
are largest near the edge of the cloud where the signal- 
to-noise ratio is low, a more effective characterization of 
the motion is through the spatial derivative of the ve- 
locity, du/dz, at the center of the cloud. The derivative 
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FIG. 3: Time evolution of the velocity gradients for temper- 
atures: T/Tc = 0.8 (■), 0.65 (•), and 0.45 (a). A typical 
predicted velocity gradient derived from the measured phase 
gradient is shown by the dashed line. The inset shows a typi- 
cal measured relative velocity profile at one point in time from 
which the gradient du/dz is derived (T/T c — 0.65). The solid 
line is the measured velocity, and the dashed line corresponds 
to the velocity predicted by the phase profile. 



contains both information about the magnitude of the 
velocity, as well as whether |1) is flowing out and |2) 
flowing in (negative slope) or vice versa. Figure [3] shows 
the time evolution of du/dz for a range of temperatures, 
T/T c . 

The first point to note in Fig. |3] is that the velocity 
oscillates in time. This oscillation is indicative of an os- 
cillation in the spin vector, i.e. of spin waves [(J, |7j. Spin 
waves exist only in the normal component, and not the 
condensate, due to the lack of condensate-condensate ex- 
change scattering. Therefore it is surprising to see signs 
of spin waves in the center of the ensemble, which is dom- 
inated by the presence of the condensate. However, it 
was shown in [111 ] that the spins of the condensate cou- 
ple to the normal component spin wave through normal 
component-condensate exchange scattering, "dragging" 
along the spins of the condensate. 

Because exchange scattering couples the condensate 
and normal component spins, the oscillatory nature of 
the velocity gradients seen in Fig. |3] implies that conden- 
sate domain formation reverses itself as well. Reversal of 
domain formation is not observed, however. There are 
several reasons this may be the case. First, the |2) state 
lifetime is approximately 150 ms for typical condensate 
densities due to dipolar relaxation. This loss shrinks the 
1 2) state population so that it does not exert a large effect 
on the |1) state and also has an apparent inward veloc- 
ity. Furthermore, the thermal component decoheres due 
to the large spatial energy inhomogeneity, as well as pos- 



sibly through spin locking with the condensate, and this 
decoherence can melt the condensate (see [10| ) . Another 
explanation for the lack of domain reversal is that the 
condensate imposes a spin stiffness to the system. It is 
energetically unfavorable for the condensate spin states 
to interpenetrate again, and phase-separated condensate 
inhibits the remainder of the spin wave cycle. 

The second point to note in Fig. [3] is that the time 
of maximal domain formation is approximately 60 ms 
(when du/dz becomes zero in the center of the sample) , 
significantly faster than domain formation time scales in 
pure condensates |4|. Spin waves are further implicated 
as the driving mechanism in the observed spin motion 
by comparing the velocities derived from the measured 
density profiles and the potential-driven velocities pre- 
dicted from the spatial phase winding. These predictions 
are made in the following way. Phase profiles such as 
in Fig. [2] are obtained using Ramsey spectroscopy [6( . A 
closely spaced tt/2 — tt/2 pulse sequence measures the dif- 
ferential phases accrued across the atom cloud due to the 
spatially varying mean-field energy and magnetic trap- 
ping potential. These phase gradients lead to velocities, 
as given by Eq. (JTJ. The dashed lines in Fig. [3] and its 
inset show typical predicted velocity profiles and velocity 
gradients derived therewith. 

It is immediately clear that there is a large discrep- 
ancy between the measured and the predicted velocities, 
which leads spin domains to form much faster than ex- 
pected. In fact, the observed velocities in the condensate 
approach the Bogoliubov speed of sound, v s — y/gn/m, 
where g = 4:Trh 2 a/m [l8|. The nature of the anomalously 
fast domain formation is compellingly shown by compar- 
ing the measured and predicted velocities (Fig. QJ . This 
comparison is parameterized by the ratio of the velocity 
gradients at the center of the cloud at the time of the peak 
measured velocity gradient. For nearly pure BEC's at the 
lowest T/T c , the velocity derived from the phase profile 
is in reasonably agreement with the actual measured spin 
velocities, but for higher temperatures the measured ve- 
locity is significantly higher than the potential-driven ve- 
locity. This difference is even more pronounced near T c , 
when the measured velocities become many times greater 
than the potential-driven motion. Lower temperature 
systems agree more closely with the value expected from 
the phase gradients, showing a lessened effect of conden- 
sate spin exchange with the smaller normal component. 
The interpretation of the velocity enhancement shown in 
Fig.|4]is the entrainment of condensate domain formation 
by spin waves in the normal component. Spin waves ac- 
celerate the condensate to velocities up to several orders 
of magnitude larger than those driven by inhomogeneities 
in the potential, and condensate-normal component ex- 
change scattering couples the condensate spin to normal 
component spin waves. 

These results agree qualitatively with the Ginzberg- 
Landau coupled-fluid theory presented in [9J . At T c they 
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FIG. 4: Enhancement of the expected relative velocity, that 
is the ratio of the measured relative velocity [Eq. J2J] to the 
relative velocity predicted by the spatially varying phase of 
the superposition [Eq. fl])]. Values significantly greater than 
1 represent a strong coupling of the condensate to normal 
component spin waves. Spin waves are suppressed at low 
temperatures when the normal component-normal component 
exchange scattering rate becomes too low, and the superfluid 
velocity approaches that predicted by the phase profiles. 



predict a sharp increase in the spin stiffness, quantified 
by the coupled-fluid stiffness parameter c s , which governs 
the spin wave dispersion relation Huj — c s k 2 for spin wave 
frequency w and wavevector k. Below T c , c s is predicted 
to be a sum of the normal component spin stiffness Cq 
and a term roughly proportional to the condensate order 
parameter. In the strong coupling limit, where the axial 
trap frequency is much smaller than the exchange colli- 
sion frequency, w eX ch = gn/U, cq is given by the exchange 
energy according to [7J 
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The explicit temperature dependence of cq is canceled by 
that of k (proportional to the size of the atom cloud), and 
a constant spin wave frequency is observed as a function 
of temperature for fixed density [fj|. However, as soon as 
there is a significant condensate component present, the 
spin wave frequency increases several fold. It is difficult 
to make quantitative comparisons with theory below T c 
for several reasons. First, as noted previously, the con- 
densate inhibits spin wave oscillations after a half cycle, 
making it difficult to determine the oscillation frequency. 
More importantly, the observed spin waves are strongly 
driven in the nonlinear regime and exhibit amplitude de- 
pendence [a 0|. Above T c it is possible to modify the 
excitation strength by changing the magnetic contribu- 
tion to the inhomogeneous potential that drives the spin 
waves and to extrapolate the spin wave frequency to the 
low-amplitude limit. Below T c , however, the presence of 
the large inhomogeneous mean-field potential of the con- 
densate makes it impossible to flatten the curvature of 
the driving field through similar means. 

In conclusion, we have measured the relative veloci- 



ties of spin components in a partially condensed spinor 
system. By comparing the measured velocities with 
potential-driven velocities predicted from phase winding, 
we rule out potential-driven superfluid flow as a mech- 
anism for domain formation in the condensate and find 
coupling between the condensate and normal component 
spin waves to be the driving mechanism for energy trans- 
port. These results agree qualitatively with the theoret- 
ical predictions. Additionally, while spin wave behavior 
drives the condensate, the spin stiffness of the phase- 
separated condensate inhibits further evolution of spin 
waves. This complex interdependent behavior highlights 
the richness of a partially condensed spinor Bose gas. 

The author thanks Ryan Thomas for experimental as- 
sistance. This work was supported by NSERC and CFI. 
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